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^ Abstract 

^— > 

o 

The dynamic surface critical behavior of macroscopic systems whose dynamic 
bulk critical behavior is described by model B is investigated. The semi- 
infinite extensions of bulk model B introduced in a previous treatment [Phys. 

CO 

f"^**) Rev. B 49, 2846 (1994)] and called models Ba and Bb are studied in detail 

o 



by means of field-theoretic RG methods. The distinctive feature of these 
models is the presence or absence of relevant nonconservative surface terms. 
The earlier results on the structure of the required counterterms, the flow 

i 

equations, and the representation of the surface critical exponents of dynamic 
q quantities at the ordinary and special phase transitions in terms of static 

bulk and surface exponents are corroborated by means of an explicit two-loop 
calculation for 4— e dimensions. For parameter values corresponding to a given 
static surface universality class, models Ba and Bb represent distinct dynamic 
surface universality classes. Differences between these manifest themselves in 
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the behavior of scaling functions for dynamic surface susceptibilities. RG- 
improved perturbation theory is used to compute some of these susceptibilities 
to one-loop order. 

I. INTRODUCTION 

The past two decades have witnessed enormous activity and progress in the field of 
critical phenomena at boundaries, such as surfaces, interfaces, or walls [1-4]. This has 
led to the development of a fairly satisfactory theory for static phenomena of this kind in 
thermal equilibrium systems. As far as quantitative predictions are concerned, the theory 
has not quite yet reached the degree of accuracy of the theory of bulk critical phenomena. 
However, practically all the sophisticated quantitative techniques of the modern theory of 
critical phenomena have been extended to systems with boundaries. Thus there are no 
fundamental reasons that would prevent one from achieving the same numerical precision as 
in the bulk case. The quality of quantitative results that can be attained through theoretical 
analyses is mainly determined by how much effort one is willing to spend [5-8]. 

An equally important and impressive development has taken place on the experimental 
side. During the past decade powerful surface-sensitive techniques have been developed and 
refined [3]. As a consequence, accurate experimental studies of static critical phenomena at 
surfaces and their quantitative comparison with theoretical predictions have become possible. 

The situation is much less satisfactory if one considers dynamic critical phenomena at 
surfaces. Despite the high level of sophistication of the available experimental techniques, 
the quantitative study of such dynamic surface critical behavior remains a major challenge. 
On the theoretical side, the pertinent knowledge accumulated is rather limited: only a few 
semi-infinite extensions of known bulk models [18] have been examined in some detail [9-17], 
and reliable quantitative results are scarce. Furthermore, from a more fundamental point of 
view one would like the theory to achieve a classification of the possible types of dynamic 
critical behavior at surfaces into corresponding dynamic surface universality classes and 



2 



their representation in terms of appropriate prototype continuum models. This analog of 
the program expounded in Ref. [18] for dynamic bulk critical behavior has only just been 
started. 

In a previous paper [17], hereafter referred to as I, the dynamic surface critical behavior 
of semi-infinite systems with the following properties was studied: 

(i) Their dynamic bulk critical behavior is represented by model B of Ref. [18]. 

(ii) Their static surface critical behavior is described by the semi-infinite <f> 4 model of 
Ref. [2]. 

Because of (i), the order parameter <f> must be a conserved density away from the surface, 
but it need not have this property in the vicinity of the surface. Such a local violation of 
the continuity equation was found to correspond to a relevant surface perturbation whose 
strength could be parametrized by a surface variable cq. The corresponding extensions of 
model B to the d-dimensional half-space z > for the cases cq > with, and Co = without, 
nonconservative surface terms were called models Ba and Bb, respectively. Both represent 
distinct dynamic surface universality classes. Hence each one of the static universality classes 
in question — namely, those of the ordinary, special, and extraordinary transitions [2] 
splits up into two distinct dynamic ones. 

In I, models Ba and Bb were investigated using field-theoretic renormalization group 
(RG) methods in 4 — e dimensions. Invoking a combination of arguments, the general 
structure of the necessary counterterms was derived. This led to the assertion that the 
required renormalization functions could all be expressed in terms of known ones. In turn, 
this implied that the dynamic (bulk and surface) critical indices can all be written in terms 
of known static (bulk and surface) exponents. 

The aim of the present paper is twofold: to verify the findings of I by explicitly carrying 
through the field-theoretic RG program to two-loop order, and to compute quantities that 
discriminate between the dynamic surface universality classes represented by models Ba and 
Bb- Obvious candidates for such quantities are the universal scaling functions associated 
with dynamic surface susceptibilities. Using RG-improved perturbation theory, we have 

3 



calculated such scaling functions to one-loop order for the case of the special transition. 
Owing to the complicated form of the free dynamic propagators, these calculations are 
rather demanding. This has kept us from computing more scaling functions to the same 
order of RG-improved perturbation theory. 

The remainder of the paper is organized a follows. In Sec. II we briefly recapitulate the 
definition of models B^ and Bg and their functional-integral representation. In Sec. Ill we 
describe details of the calculations needed to verify the renormalization of the theory to two- 
loop order. In Sec. IV our one-loop results for susceptibility scaling functions are presented. 
This includes a discussion of their crossover from critical to hydrodynamic behavior. Section 
V contains a brief discussion of our results and conclusions. In Appendices A-F various 
calculational details are explicated. 

II. DEFINITION AND FUNCTIONAL-INTEGRAL REPRESENTATION OF 

MODELS B A AND B b 

We consider a time-dependent n-component order-parameter field 0(x, t) = (0 Q (x, t)) on 
the half-space V = {x = (xy, z) | X| £ iR d_1 , z > 0} bounded by the surface dV = {(xy, z = 
0) | X|| £ IR^ 1 }. As expounded in I, the models we are concerned with are defined by the 
Langevin equation 



is the Hamiltonian of the semi-infinite (j) 4 model. The operator D is defined through its 
action 





Here Aq is a constant, and 







on functions / over V satisfying the boundary condition 
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d n f = £of (4) 

on dV, where the derivative d n is along the inner normal. This Robin boundary condition 
[19] guarantees that the self-adjoint extension of D exists; it is satisfied by the functions 

^(x, t) = [(-A + r + \ u <f> 2 ) <f>] (x, t) (5) 

produced by the functional derivative SH/5(f) in (1). 
The Gaussian noise ( has mean zero and variance 

(C Q (x,t)C /3 (x / ,0) = 2A 5 Q/3 D^8(t-t') , (6) 

where D X x' is the integral kernel associated with D. Writing it in the formally symmetric 
form 

D xx , = W(x - x') V' + £ S(z) 5(x - x') , (7) 

one sees that the variable cq measures indeed the strength of the nonconservative surface 
terms. 

The above equations with cq > and cq = define models Ba and Bb, respectively. 

In order to apply renormalized perturbation theory, it is convenient to transform to the 
equivalent functional-integral representation. According to I, this gives the dynamic action 
[17] 

J = j dt [j v [^-A o (A0)(^-0)] 

+ A | [(d n - cq)$\ [H^-4>- 8(0)(d n - c o )0] + (A0)(d„ - c )^ , (8) 

in which <p is the usual auxiliary field [20]. The singularity oc 8(0) is a boundary term that 
arises from the surface potential oc 5(z) in (2) upon integration by parts. To avoid it, one 
can replace the latter ^-function by an appropriate regularized one, such as 

8 B {z)=Be- Bz (9) 
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with arbitrarily large but finite B. Then 5(0) gets replaced by B. Since the boundary 
conditions found in I, 

d n <f> = c Q <f> , dnU^ = cqH^ (10a) 

d n A(j) = c A(j) , d n (j) = c 4> , (10b) 

imply that the contributions from the surface integral in (8) vanish, we do not have to worry 
any further about the singular term oc 5(0) in the action J . 



III. RENORMALIZATION 

As in I, we consider the general response and correlation functions 

IN N M M \ C 

n(-AoA^) n<^n(- wi 8 ) n 



s I 1 



, j=l k=l 1=1 m=l 



fill 



as well as the related functions 

G (N,N;M,M)^ t) 

IN N M M \C 

= (n*n^n*i. n*i.) ■ (12) 

\ j = l k=l 1=1 m=l I 

Here the superscript C indicates cumulant averages. The variables x, r, t represent the 
N + N points off the surface, the M + M parallel coordinates of the surface points, and the 
set of all time arguments, respectively. 

According to the arguments given in I, to renormalize these functions in 4 — e dimensions, 
it should be sufficient to supplement the usual static reparametrizations [2] 

</> = Z 1 ' 2 ^ , T = ^Z t t + n , (13a) 



u = !Lz u u, s d =(47i)- d ^\ (13b) 

Sd 



</>\ 8 = (Z <t> Z 1 ) 1 / 2 (<l>\ 8 ) R , c = iiZ c c + c sp , (13c) 
and the well-known dynamic ones 

4> = Z- 1/2 4> R , \ = /i- 4 Z^\;, (14a) 

by 

4>\s = Z^ 1/2 (4>\ S ) R , c Q = fie , (15a) 

(A$) S = (Z 1 /Z^ 2 [(A$) S ] R . (15b) 
Thus the renormalized functions [21] 

G (N,N;M,M) = z ^^ z -f q(N,N;M ,M) (1 g } 

and 

iiTr »r iTr i,r\ N+M+N+M M+M 

W (N,N-MM) = z ~ t z --^~ W (N,N;M,M) 



should be ultraviolet-finite. Utilizing dimensional regularization, we will verify this to the 
order of two-loops for the functions Q^' 1,1,0 ^ an d W R 0,1 ' 1,0 \ To this end we shall first compute 
these functions for r = c = 0. In a second step we shall then check the cancellation of c- 
dependent poles by expanding about the special point. For notational simplicity we will set 
the number of components n = 1 henceforth. 



A. Graphical notation 

The Feynman graphs of the functions we are interested in involve the free propagators 

G(x,t;x',O = (0(x,t)^x',t'))/ (18) 

and 

C(x,t;x',0 = <<^M)0(x',t')}/ . (19) 



It is favorable to work in a pzoj representation. For their Fourier transforms, defined by 



G(p;z,z ;u) = 



jd^x^j dt e - i[p - (x " - x i|)^(*-*')] G(x, t ; x' , t') 



(20) 



we use the graphical notation 



O- 

z 



p, CO 



(21) 



and 



_ p, co _ 

° °< = C(p; 2, 



(22) 



In order to indicate whether external points, such as the above ones labelled z and z' , are 
off or on the surface, we use open circles (as shown) for points with z > and crossed circles 
(as in Fig. 1) for surface points. 

The explicit expression for the response propagator (21) may be gleaned from I; since it 
is rather lengthy, even in the special case r = c = 0, we have relegated it to Appendix A 
[see Eq. (Al)]. 

The correlation propagator C can be expressed in terms of G, exploiting the fluctuation- 
dissipation theorem. One has 



o O 2A T 

z z' = Im 



p, CO 

O ~* *-O r 

z z 



LO 



Im[x [0] (p;^,^;^; 



(23) 



where the full triangle stands for the negative Lapacian p 2 — d 2 z , and i s a zero-loop 
susceptibility. 

In conformity with these conventions, we represent the four-point vertex of the action 
(8) by 



(24) 
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As a further graphical element, we shall need the vertex representing insertions of the 
operator — Aq / dt j dv <f>A<f> to which cq couples. We use the symbol 



(25) 



B. Calculation for c = 

In Fig. 1 the graphs of G^ 0,1 ' 1 ' ) are depicted up to two-loop order. The corresponding 
graphs of VL^ ' 1 ' 1 ' ) /Xq are obtained from these by putting a full triangle on one external leg. 

To compute the latter graphs, it is convenient to proceed as follows. We amputate the free 
propagator lines attached to the external points. The resulting amputated graphs must then 
be considered as distributions with respect to their dependence on z and z' . Regularized by 
analytic continuation in e, these distributions are well-defined and can be Laurent expanded 
in e [22]. For simple powers such as z~ 2+€ , the Laurent expansion is well-known and can be 
looked up in Ref. [22]. However, in general one must study the action of the distributions 
on test functions to determine their e expansion. To illustrate the procedure, let us consider 
the integral 

/•no 

dzV e (z)f(z), (26) 



o 

where T> t {z) is a distribution and f(z) a test function. As a test function, f(z) can be 
written in terms of its Laplace transform f(a) as 

C-Hoo 

f(z)= Jdae—f(a), (27) 

C — ioc 

with arbitrary C > 0. Hence it is sufficient to know how T> £ acts on the class of functions e~ az 
with Re a > 0. Upon computing the corresponding (dimensionally regularized) integrals, 
one arrives at Laurent series of the form 

dzV e (z)e—= £ ^-fi, (28) 

j--n 



where uq is the order of the highest pole. In the next step one has to find distributions 
T)W (z) that solve 

roc 

/ dzV (i \z)e- az = d {j) (a) . (29) 
Jo 

This yields the desired expansion 

V e (z) = £ ^-p- . (30) 

j--n a 

In Appendix A the e expansion of the amputated two-point graphs resulting from the 
graphs of Fig. 1 are given to the required order in e. From these results the corresponding 
e expansion of each W and G-function with N + N + M + M = 2 can be obtained in 
a straightforward fashion by applying the distributions to the respective products of free 
propagators that the amputated external legs represent. The explicit expressions one obtains 
for the graphs of G^ ' 1 ' 1,0 ' and I4A ' 1 ' 1 ' ) are given in Appendix B. 

To determine the renormalized functions, we substitute these results into (16) and (17) 
and use the reparametrizations (13a)-(15b) together with the known two-loop results [2] 

Z<p= l--^u 2 + 0(u 3 ) , (31) 
Z u = l + ^u + 0(u 2 ) , (32) 



Z 1 = l + -u + 

e 



2 1 
7 ~ 7 



u 2 + 0(u 3 ) (33) 



for n = 1. Proceeding in this fashion, we have explicitly verified to two-loop order the 
cancellation of all dimensional poles in G^' 1,1 ' ' and IT^ ' 1 ' 1 ' ^. Details of these lengthy 
calculations are given in Appendix C. Their final results read 

Gf hl ' 0) (p, z, u; A, c, u) = — 9 -r— - ) — — \2k + e~ k - z - 2k_ e~ k+z 

K v ; \(k% - k 2 _)(2k + k_ + c(k + + I + 

- \u[S + e- k+z -£_e~ k - z } + u 2 0(e°) + 0(u 3 )} (34) 

and 

<'«(p, A, c, u) = 2k+k+ l {k++k) {[(k- + ^ + (*+ + c)e-- 
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+ \u(\£ E (k_ + c) - \T + ]e~ k+z + [C E (k + + £)- \T_ 
+ « 2 O(e ) + O(« 3 )| . 



(35) 



Here C E is Euler's number, and the quantities £±, T±, and k± are defined by Eqs. (B3) and 
(C5), by (BIO) and (Cll), and by (C3) of Appendices B and C, respectively. 

C. Expansion about the special point 

Until now our calculations were restricted to the case of critical surface enhancement, c = 
0. We now wish to check whether the chosen reparametrization also absorb the additional 
poles that appear for 0. In order to bypass the complicated computation of c-dependent 
graphs, we expand about the multicritical theory Since the critical bare surface enhancement 
c sp vanishes in our perturbative approach based on dimensional regularization, this amounts 
to an expansion in powers of cq for the bare theory. For simplicity, we will restrict ourselves 
to the case of model B E in this calculation, setting c = 0. 

From the action (8) we see that the operation —d/dcQ produces an insertion of the (Fourier 
transformed) operator — \q[(4>A4>) s ]p ,n a ^ zero parallel momentum P and zero frequency Q. 
Let us denote the analogs of the functions q( n ^ n 'MM) w ith / insertions of — Ao[(0A0) s ]p^ 
by G( JV ' JV ^ M ' M ; / ) . Then the expansion of the functions G (or W) in powers of cq can be 
written as 

q(N,N\M,M) _ q(N,N;M,M) 

_ Co £(N,N;M,M;l) +0(c 2) ? (36) 

where the functions on the right-hand side are taken at cq = 0, and it is understood that 
the inserted operator has P = and $7 = 0. 

A word of caution must be added here: We wish to use this expansion in the massless case 
t = 0, where it becomes problematic because of infrared singularities [23]. Strictly speaking, 
we should take the limit P — >• of the massless functions G( jV ' jV ; M ' M ; / ) only after the infrared 
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singularities have been resummed via the RG [24]. For the graphs of & N ' N ' M < M 'V we are 
going to consider, setting P = causes no problems as their ultraviolet poles can be safely 
identified. 

According to the reparametrizations (13a)-(15b), the renormalized G-functions with 
insertions should be given by 

, - , . - , . T , N+M-N-M M , ~ ~ 

G (N,N;M,M;D = Z ^~ 2 Z -~ Z I G (N,N;M,M;I) ( g ?) 

We will verify to the order of two loops that all ultraviolet poles cancel in G^' 1 ' 1 ' ' 1 ^, if we 
use the previously given two-loop expressions for the renormalization factors together with 
the known n = 1 result [24,2] 




The graphs of G^°' 1;1,0;1 ' contributing to this order are shown in Fig. 2. Appendix D 
contains the e expansions of the required distributions and of the graphs themselves. Upon 
making the transition to renormalized quantities and summing the contributions, one arrives 
at 

G^ 1 ;W) (P, *, A, c = 0,u) = + u 0(6°) + , 2 0( e °) + 0(u 3 ) . (39) 

IV. SCALING FUNCTIONS 

According to (15a) cq is not renormalized. This implies the trivial flow behavior 

m = y (40) 

under RG transformations \x —¥ fil. On physical grounds the initial value c must be non- 
negative. Hence there are two fixed-point values of physical interest: c* A = oo and c* B = 0. 
In the enlarged parameter space of the dynamic theory, each static fixed point (describing 
the ordinary, special, or extraordinary transition) unfolds into a line c > between pairs of 
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dynamic fixed points with these two values of c. The fixed points corresponding to models 
Ba and Bb are infrared-stable and infrared-unstable in the c-direction, respectively (see the 
flow diagram depicted in Fig. 5 of I). 

The members of each of these pairs of fixed points represent, distinct dynamic surface 
universality classes with the same static critical behavior. Since the dynamic critical expo- 
nents can be expressed in terms of static ones, we must look at other universal properties 
to see the difference between these universality classes — and hence between the dynamic 
surface critical behavior of model Ba and of model Bb- 

As a simple illustrative quantity we will examine the surface susceptibility 



Xn(r,<) 



5/ii,o(0,0) 



= W (p > 1 >°> 1) (r,t;O,0) (41) 



for c = and r > 0. Here hi^(r,t) is a position and time-dependent surface magnetic field. 

From the general results for the VF-functions given in Eqs. (3.67) and (3.69a)-(3.70) of I 
we can infer the scaling form of Xn( r i ^) f° r r > and c = 0. (As usual, these scaling forms 
are obtained in our field-theoretic approach as fixed-point solutions to the RG equations for 
the renormalized functions, but they describe the critical behavior of the bare functions as 
well. We therefore drop the subscript R in the rest of this section.) 

Upon Fourier transformation we arrive at the scaling form 

Xn(p,6j, r,c) «Cir-T"*(p,w,c) (42) 

with 

V = C 2 pT~ v , w = C 3 cot- v( , c = C 4 ct- v . (43) 

Here 7^ = v(\ — rf^) is the usual surface susceptibility exponent of the special transition, 
C = 4 — rj denotes the dynamic bulk exponent, and the C; are nonuniversal constants. 

As is borne out by our perturbative results (see Appendix E and below), the left-hand 
side of (42) exists for c = c* A and c = c* B . Hence the corresponding limiting values 

13 



*a,b(p,w) = *(p,w,£* ) 



(44) 



of the scaling function ^ for c — Y c* A = oo and c — Y c* A = should exist. Utilizing RG- 
improved perturbation theory, we have computed the scaling functions and to one- 
loop order. Some details of the calculation can be found in Appendix E. The results are 

2(k+ + k_ 



*MP,w) 



+ 



(k + + k_) 2 + l 12k+k_(k+ + k_)(k+ - k_) 2 [(k+ + k_) 2 + l] 2 
■\A + B + + B_ + 4k+k_(k* - k 2 f + (1 - C E )(k + - k_) 2 



(k+ + k_) 2 (k+ + 10k + k + k 2 ) - 2(k + + k_) 2 + 1 



+ 0(< 



(45) 



and 



*b(p,w) 



(k + + k_) 2 -l 

2k+k_(k+ + k_) ' 48k^ k 3 _ (k+ - k_) 2 (k+ + k_) 3 



+ 



• \ k+k_ A + ki B + + kl B_ - 4ki kl (k 



2 i 2 /i 2 i 2 \2 



+ (l-C £ )(k + -k_ 



4 /i,2 



(k+ + k_) 4 (k 2 - k + k_ + kl) 



2(k + + k_) 2 (k 2 + + k + k_ + k 2 _) + (k 2 + + 3k + k_ + k 2 _) } + 0{e 2 ) , (46) 



with 



A = 4k+k_ 



'k+ + k_ - 2 
k+ + k_ + 2 



In 



Vk + + k + 2 + Vk + + k_ - 2 
Vk+ + k_ + 2 - Vk+ + k_ - 2 



(k+ + k 



+ 2(k+ + k_) - (k+ - k_) 2 (k+ + k_) 4 + 2(2k^k!. + 2k 2 + k 3 _ -k 5 + - k 4 k_ - k+k 4 



+ (2 - C 1 



k 6 + 2klk + - kiki - 4k 3 _ki - k 2 k 4 + 2k_ki + ki 



(47) 



and 

B± = 

where 



-4k ± k T (k 2 - l) 1 / 2 (k ± + k T ) (k| - k 2 - l) 2 In | 



' y/k^TT+ Vk~^T ' 
,v / I%TT-v / k^T. 



(4? 



k± = 



\ 



p 2 + - ± \ - + iw 



(49) 
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A well-known general property of our model is that dynamic susceptibilities .) in 
the limit uj — > go over into their static counterparts % stat (.). Hence 

lim x n (p, u, r, c, c) = xff(p, r, c) . (50) 

This implies that as w — >■ with all other arguments fixed, the scaling function ^(p,w,c) 
— and hence both ^ as well ^>b — must approach a limiting function that is independent 
of c and equal to the corresponding scaling function of x| t 1 at . Setting again c = 0, we write 

^(p^jKdr-^^'lp) . (51) 

Then 

£(p,w = 0,c) = ^ A (p,w = 0) = ^ B (p,w = 0) 

= * stat (p) . (52) 

As an explicit check the reader may verify that for w = the e expansions (45) and (46) of 
^>a and ^ b both reduce to the same result 



Vjjstat/ \ = 1 — 1— 1 

iPJ vT^TT + 3 



2 — Ce 1 — Ce 



2 V / P T TT 4(p2 + l) 3 / 2 
-^ T ln (p+^ + T )]+°( e2 ) (53) 

of the static theory. 

On the other hand, the functions ^>a and ^>b differ in general whenever w^fl. Inspection 
of their w-dependence reveals, in particular, that 

Re [^ A (p = 0,w^0)] - vj) stat (0) - w 3/2 , (54) 

but 

Re [\& B (p = 0,w -►())] - vj) stat (0) - w 1/2 . (55) 

To elucidate further these differences between the corresponding dynamic surface uni- 
versality classes of models Ba and Bb, we consider the behavior of \n as a function of t. 
Fourier transformation of (42) yields 
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Xii(p,t,r,2)^C 1 C 3 - 1 r^-^*(p,t,£) 

= C 1 C 3 - 1 t-( c - 1+ ^H(p,t,c) (56) 

with t = C^tr 1 "* , where \&(.,t, .) is the inverse Fourier transform of ^(.,w, .). Depending 
on whether c is small or large, we may set c to the fixed-point values c* B = and c* B = oo. 
By analogy with before, we denote the resulting scaling functions by E b , ^a, and Ea, 
respectively. 

We could perform the Fourier inversion only numerically. The numerical results for 
x 5 J 4(p = 0,t) and ^ B (p = 1 t) one obtains with e set equal to one are plotted in Fig. 3. Since 
at criticality, the p = susceptibility Xn has a long-time tail ~ t v \\ , the limiting 
values 5*4 = 5^(0,0) and E B = Eg (0,0) both exists, which implies the small-t behavior 

^,B(p = 0,t->0) «H° iB r (f - 1+ ^ /f . (57) 

Hence the difference between the dynamic surface universality classes manifest itself here 
only through different values of the universal amplitudes E A and E B . 

The asymptotic forms of ^a,b(0,^) in the large-X. limit reflect the time-dependence of 
Xn (p = 0, t, t > 0, c) for c > and c = 0. Since the order parameter is a conserved density 
away from the surface, one anticipates a hydrodynamic long-time tail even in the case of 
model Ba- The asymptotic small-frequency dependencies given in (54) and (55) suggest 
distinct algebraic decays ~ f -5 / 2 and ~ f -3 / 2 for models Ba and B B , respectively. 

In Appendix F the large-t behavior of the zero-loop scaling functions (p, t) is worked 
out in detail. We find the asymptotic forms 



(P ' tj " 2v ^(l + p 2 ) 2 

and 

yTTW p2(1+p2)t 

-B (P.*) = 75 IT/i 



, 3(5p 2 + 5p 6 - 2) 1 / 1 
P ~ ,« . 7^77-^7772 7 + 77 



4(l + p 2 ) 2 (l + 2p 2 ) 2 t Vt 



, 2 + p 2 1 /l 
P + 4(l + 2p 2 ) 2 T + lF 



(5? 



(59) 



For p > 0, the results exhibit the expected exponential decay with the same relaxation rate 
p 2 (l + p 2 ). In the hydrodynamic limit p — > 0, the relaxation rate vanishes, and we get the 
algebraic decay laws 
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^ A (p = 0,t^oo) - 1T 5/2 (60) 

and 

^ B (p = 0,t^oo) - 1T 3/2 . (61) 

Since the above large-t behavior is of a noncritical, purely hydrodynamic, origin, its form 
should already be correctly obtained at zero-loop order. Terms beyond this order will affect 
the (suppressed) universal amplitudes involved in (60) and (61), but should neither change 
the exponents in these equations nor the exponential decay on a time scale oc p -2 predicted in 
(58) and (59). Our one-loop results shown in Fig. 3 are in conformity with these expectations. 

V. CONCLUDING REMARKS 

Let us conclude by summarizing the principal results of the present work. 

Using field-theoretic RG methods, we investigated the dynamic surface critical behavior 
of the semi-infinite models B^ and Bg introduced in I and Ref. [16]. These models represent 
the dynamic surface universality classes of surface-bounded macroscopic systems whose dy- 
namic bulk critical behavior is described by the familiar model B of Ref. [18]). They differ 
by the presence or absence of nonconservative surface terms. 

Our first goal was to check whether the dynamic field theories of these models can 
indeed be renormalized in the manner asserted in I. To this end, we presented a variety of 
two-loop calculations through which we confirmed the claim of I that all renormalization 
factors involved are given by the known renormalization factors of the static theory. The 
resulting RG equation then lead to the conclusion that the critical exponents of the dynamic 
theory are expressible in terms of static (bulk and surface) critical indices. Since models Ba 
and Bb have the same thermodynamic equilibrium state, the dynamic critical exponents 
of both models are the same and hence do not differentiate between the respective two 
distinct dynamic surface universality classes that are associated with the same static surface 
universality class. 
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Our second goal was to investigate a characteristic dynamic quantity that is sensitive 
to the differences between models and Bb, and hence distinguishes between the cor- 
responding dynamic surface universality classes. Using RG-improved perturbation theory, 
we computed the dynamic surface susceptibility xn(p,w) for the static universality class 
of the special transition (c = 0) and the associated scaling functions a an d ^b- The e 
expansions of the latter functions are given in (45) and (46). Their analogs ^^(C^t) and 
\&B(0,t) for Xn(p = 0,t) are depicted in Fig. 3. They are clearly different; in particular, 
their asymptotic behavior for large values of the scaling argument t ~ tr 1 "^ is markedly 
different. 

As our calculations show, detailed studies of dynamic critical behavior at surfaces by 
means of field-theoretic RG methods and the e expansion, including the calculation of scaling 
functions, are feasible, albeit technically rather demanding. Unfortunately, we are at present 
not aware of experimental studies of dynamic surface critical behavior of systems belonging 
to the bulk universality class of model B. Once such studies become available, the measured 
observables could be computed by the methods described above. 
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APPENDIX A: e EXPANSION OF DISTRIBUTIONS ASSOCIATED WITH 
AMPUTATED TWO-POINT GRAPHS 

In the following, r and c are set to their (multi-) critical values r = c = 0. The response 
propagator G then becomes 

O ^ o 1 f 1 2k + k_+c {k + - k_) {z+zl) 



, & — G 
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2k + k_ + cq (k + + k. 



k_[2k + k_ + Cq (k + + 



with 



1/2 



We also need the well-known static propagator 



e -p(z+z') _|_ e -p\z-z'\ 



(Al) 



(A2) 



z Jstat 2p 

The required one-loop graph is proportional to its static counterpart. We have 
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(A6) 



Here Ce denotes Euler's number, and z+ m+€ is a distribution whose definition for e > can 
be found in Ref. [22] or in the appendix of Ref. [2]. 
In a similar fashion we get 
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The graph that contains the static one-loop graph (A4) twice is not simply proportional 
to a static one. The e expansion of this distribution can be determined by computing its 
action on the test function e~ az+t3z ' (cf. Sec. III. A). This yields 




— ^— |- [5(z) 8(z) + 8(z) 8(z ) + 2K+K _ + + K _ } S (z) 6 (z ) 
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+ 0(6° 
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The Laurent part of the e expansion of the remaining two-loop graph may be conveniently 
determined by expanding in powers of u. It is easy to see that the ^-dependent part of the 
graph is regular in e. (Use power counting and the form of the divergent subgraphs.) Hence 
the poles result from the co = term, which is proportional to the corresponding graph of 
the static theory. It follows that 
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APPENDIX B: GRAPHS OF G(°' 1;1 ' ' AND W^ ' 1 ' 1 ® FOR r = c = 



Utilizing the results of Appendix A, the e expansion of the desired graphs of G^ ' 1 ' 1 ' ) 
and W^ 1 '^ for r = c = can be derived in a straightforward fashion. We first list our 
results for the graphs of G^ ' 1 ' 1 ' ). To one-loop order we have the contributions 
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where Ei is the exponential integral function. 
Our results for the two-loop graphs are 
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and 
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For the graphs of VI^ 0,1 ' 1 ' ) we find 
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APPENDIX C: REPARAMETRIZED GRAPHS OF G 



(0,1;1,0) 
R 



AND W 



(0,l;l,0) 



In the foregoing appendix the graphs of the bare functions G^' 1,1 ' ' and H 7 ^ 0,1 ' 1 '^ were 
obtained. Here we express these results in terms of the renormalized variables u and A, 
employing the reparametrizations (13b) and (14a). For notational simplicity we set \i = 1. 



We need the n = 1 relations 



An = A 
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Since we restrict ourselves to two-loop order, we must retain only terms to order w 2 . 
Hence for the two-loop graphs of Appendix B, the reparametrization amounts to the sub- 
stitutions uqSj — y u, An — > A, and Cn -4- c. Accordingly we list here only the results for 
the zero-loop and one-loop graphs. To indicate that the graphs are expressed in terms of 
renormalized parameters, we use the notation [.]R ep - 

Upon introducing the quantities 
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our results for the graphs of G^' 1 ' 1 ' ^ may be written as 
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Our results for the graphs of W^ ' 1,1 ' ' are 
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with 

F± = F ±\ K± ^k ± ,co^c ■ ( C11 ) 

APPENDIX D: AMPUTATED GRAPHS WITH AN INSERTION OF(<f>A4>) s AND 

GRAPHS OF G(°' 1;1 ' 0;1 ) 

The results listed below were obtained by following our previous strategy: We first 
computed the Laurent expansion of the distributions associated with amputated graphs of 
(2(o,i;i,0;i) ; re l a ting them to expressions known from the static theory whenever possible. 
Subsequently we utilized these results together with those of Appendices A-C to compute 
the Laurent expansion of the graphs of G^ 0,1 ' 1 ' ' 1 ) to the required order in e. 

The graphs of the static theory appearing in the formulae below involve an insertion 
of the surface operator [§(</> 2 )s]p, the static analog of — Ao[((/)A<^) s ]p i n- As discussed in the 
main text (cf. Sec. III.C), we have set P = and Q = 0. 

It should also be recalled that we restrict ourselves to the case of model Bg, setting 
Co = in this appendix. Accordingly the functions F± and G± appearing in the results 
listed below must be interpreted as the expressions to which (BIO) and (B7) reduce for this 
particular choice of cq. 

Our results for the amputated graphs are 
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We now turn to the graphs of GK°' 1;1 ' 0;1 ). In addition to the above results and those of 
Appendix A, we need the graph 
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As an immediate consequence we have 
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For the one and twodoop graphs we find the results 
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\q(k\ — k 2 _)k + k 

Upon insertion of the relations (CI) for uq and Aq one obtains the reparametrized graphs. 
The explicit expressions resulting to order u 2 for the zero and onedoop graphs read 
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APPENDIX E: CALCULATION OF SUSCEPTIBILITY SCALING FUNCTIONS 

In this appendix we give details of the calculation of the scaling functions ^a(p, w ) an d 
x I'b(p,w) whose e expansions are given in (45) and (46), respectively. In addition we explain 
how the asymptotic long-time behavior quoted in (60) and (61) can be deduced. Unlike the 
other appendices, we now have r > rather than r = 0. Accordingly the quantities k± now 
stand no longer for the tq = expressions (A2) but for 
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From the zero-loop results of Refs. [16] and [17] one obtains 
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The one-loop contributions to xn involve the vertex graph 
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a result known from the static theory and independent of cq. Here Ki_ f /2{x) is a standard 
Bessel function. 

Utilizing this result and performing the necessary integrations, one arrives at 
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where F e is the hypergeometric function 
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We use the expansion 



1 + ai 



€ _V^ ln l±vg + Q(€2) 



(E8) 



(1 - x) 2 (1 - x) 2 ~ 1 

With the aid of the the reparametrizations (13a)-(15b) and the perturbative results for the 
Z-factors quoted in (31), the renormalized susceptibilities Xii,r(c = 0) an d Xii,r(c = °°) 
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then follow in a straightforward fashion. Evaluating these at u = u* = e/3 + 0(e), the value 
of the infrared-stable fixed point, leads to the results for an d presented in (46) and 
(45). 



APPENDIX F: LONG-TIME BEHAVIOR OF SUSCEPTIBILITY SCALING 

FUNCTIONS 



We wish to determine the asymptotic behavior of the zero-loop scaling functions 



*%(P>t) = r dw*° (p,w)e- 



-«Wt 



(Fl) 



27T j-oc 

as t — y oo. The functions ^ A B (p,w) have a branch cut for Imw < — p 2 (l — p 2 ) and no pole. 
We deform the integration path as shown in Fig. 4. The contributions from ii, ^3, and ^5 
vanish. The contributions from I2 and 1 4 add up to 



*5U(p,t) = -- 



7T Jp2(l_p2) 

Expansion of Im^^ B in y = x — p 2 (l — p 2 ) yields 



dx \mi> A B (p, — ix) e xt . 



(F2) 



Im*>, -i(y + p 2 + p 4 )] = - ( y^F m ^ (P2+P4)t 

(1 + p 2 ) y 1 -' 2 



(5p 2 + 5p 6 - 2) ?/ 
2(l + 2p 2 ) 2 (l + p 2 ) 2 



+ 0(y 2 



(F3) 



and 



M>>, + p 2 + p 4 )] = - ^1 + 2p 2 



-(P 2 +P 4 )t 



(2 + p 2 ) y 



2(1 + 2p 2 



(F4) 



Upon substitution of these expansions into (F2), the integrations can be performed. This 
gives the asymptotic forms (60) and (61). 
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FIGURES 




1. Graphs of G^ - 1 ' 1 ' ). 



39 



